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1. Chaotic dynamics 

The fluctuation relation is a general symmetry prop- 
erty of mechanical systems which should hold under the 
only assumption that the system motions are chaotic: it 
reflects the time reversal symmetry. 

Time reversal symmetry means a smooth isometry / 
of phase space which anticommutes with the time evo- 
lution map St', i.e. I St — S-tl- Therefore the familiar 
operation of time reversal T or, even better, TCP would 
be always valid in fundamental models, and therefore 
the time reversal symmetry can be an issue only if one 
deals with phenomenological models in which dissipation 
is empirically introduced: as in the case of Navier-Stokes 
equations for fluids, for instance. 

The discovery, published in the paper 2] , of the fluctu- 
ation relation has led to renewed efforts towards the for- 
mulation of a theory of nonequilibrium stationary states. 
In the paper a link is attempted with the earlier proposal 
[3| for the description of the probability distribution for 
chaotic stationary states in fluids. Although the paper 
was a real breakthrough, the original argument needed to 
be made precise. A direct connection with B was estab- 
lished and called "fluctuation theorem" in [J]: where it 
was shown how the paradigm of chaotic evolution consti- 
tuted by the hyperbolic (also called "Anosov") systems 
allowed for a precise formulation of sufficient conditions 
under which the fluctuation relation held. 

The latter proof has been considered interesting be- 
cause in a sense the hyperbolic evolutions perform for 
chaotic systems the role plaid by the harmonic oscilla- 
tors for ordered systems. From these works emerged the 
interest for Physics of two new fundamental concepts, 
the chaotic hypothesis and a general mechanical notion 
for the entropy creation. 

The first is an extension of the proposals, [3j|, that iden- 
tified the probability distributions forming the ensembles 
suitable to give the statistical properties of states of tur- 
bulent fluids with the special class of distributions, well 
known and studied in the theory of dynamical systems, 
called SRB distributions. 

The extension can be formulated, [H, 0| , in the form a 
hypothesis: 

Chaotic Hypothesis: Motion on the attracting set of a 



chaotic system can be regarded as "hyperbolic". 

This hypothesis is very ambitious as it should be seen 
as an extension of the ergodic hypothesis, which it implies 
if applied to a system which is isolated and subject only 
to conservative forces. 

It allows us to state existence of time averages of me- 
chanical observables x — ► F(x), identified with func- 
tions of the phase space point x representing the micro- 
scopic state of the system, i. e. the existence of the limits 
limT^oo y Y^j=o F(S : 'x) for all initial data x chosen in 
the vicinity of the attracting set, setting aside a set of 
phase space volume. It also implies that, the limit is 
independent of x (apart from the zero volume possible 
exceptions) and therefore define a statistics, i.e. a prob- 
ability distribution \x such that 
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(1.1) 



In other words for chaotic systems motions have a well 
defined statistics, i.e. a probability distribution that al- 
lows us to define (in principle) the time averages of the 
observables. 

The nontriviality of the above statements becomes per- 
haps more clear in the case of systems subject to steady 
dissipation. Their stationary states cannot be described 
by statistics which are not singular, i.e. which attribute 
probability 1 to sets of phase space volume. 

This is a seemingly odd situation: we are interested 
in data randomly chosen with a probability distribution 
with density on phase space and, yet, they evolve with 
a statistics which is singular. Such a situation, however, 
appears to be quite clearly correct as soon as simula- 
tions are attempted in virtually any system which ex- 
hibits chaos (i.e. positive Lyapunov exponents) and is 
subject to some kind of dissipation: hence it is simply 
accounted for by the chaotic hypothesis, see for instance 
i. 

However technically the hypothesis is far more rich of 
implications. In fact the distribution fi is, in hyperbolic 
systems, identified with the SRB distribution, which is 
not well known but it should be viewed as a general- 
ization of the microcanonical distribution for isolated 
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systems. In a sense that can be made very precise it 
is the distribution that gives equal probability to cells 
into which the phase space can be imagined discretized, 
[E 0, HI ■ It allows to give a mathematical formula for the 
averages in (1.1) and even to give a precise definition of 
coarse grained cells in phase space, [g, [l(| • 

Without entering here in more technical details this 
means that we have, from the theory of the SRB distri- 
butions for hyperbolic systems, expressions for the aver- 
ages in terms of mechanical quantities. The latter can 
therefore be used to derive general relations between ob- 
servables averages even though, as it is virtually always 
the case, we cannot hope to compute their actual vlaues. 

Just as in equilibrium where we can write the aver- 
ages as integrals with respect to the Liouville distribu- 
tion on the energy surface (or the canonical one) but we 
can hardly compute them: nevertheless we can establish, 
by using the formal expressions, general relations which 
turn out to be extremely interesting precisely because of 
their generality. A celebrated example is Boltzmann's 
heat theorem, i.e. the second law of equilibrium ther- 
modynamics, as a consequence of the assumption that 
the statistics of motion (of an isolated system subject to 
conservative forces) is the microcanonical distribution. 

The chaotic hypothesis has a rather general conse- 
quence which should be seen as a generalization, at any 
distance from equilibrium, of Onsager-Machlup fluctua- 
tions theory near equilibrium, flU Ha [ri| . 

Consider a mechanical system of particles described by 
a generic equation of motion for the representative point 
x in phase space 

x = f{x) (1.2) 

where x denotes the position and velocity components. 
The equation will be a model describing a finite system 
on which external non conservative forces act. 

Therefore the equation will be non Hamiltonian and 
phase space volume will not be conserved. As a conse- 
quence the divergence 

a(x) = -J2d Xi fi(x) (1.3) 

i 

does not vanish. However in presence of time re- 
versal symmetry it will be odd under time reversal 
o~(x) = —o~(Ix). The system will be called dissipative 
if even the time average c+ = limT^oo ^ Jq o~(Stx)dt = 
J (i(dy)a(y) of a is positive. 

In a dissipative system, a + > 0, in the stationary 
state described by the statistics /z, consider the prob- 
ability that fj(Stx) ~ <p(t) for t 6 [— |r, |r] where 
t — » ip(t) is a prescribed pattern. The symbol ~ means 
that \f(S t x) — f{t)\ < s for some very small e (see 
H, EH for a quantitative form of the notion of "very 
small"). Define the dimensionless phase space contrac- 
tion o~{x) as the divergence of the equations of motion 



changed in sign. Suppose that the average (with re- 
spect to the statistics /x, i.e. the integral of o~{x) with 
respect to /x, see (1.1)) phase space contraction a + is 

positive a + > and define p — ^ f\ ^^^-dt. Sup- 
pose that fi(x), . . . , f n (x) are n observables with defined 
parity under time reversal, suppose odd for defmitcness: 
Mix) = -fi(x). Then 

Fluctuation relation Suppose that the n observables 
form a complete set in the sense that the average phase 
space contraction p is determined by the patterns followed 
by the observables f\ , . . . , /„ (for instance a is one among 
the fj or p is the sum of the averages of some of the fj). 
Then if o~+ > and the time evolution is reversible there 
exists p* > 1 and 

P T (for all j, andt £ \-\t, |r] : fj(S t x) ~ ipjjt)) _ 

P r (for all j,andt e [-±r, ±r] : fj(S t x) <fj(-t)) 

= e p CT+ T+o(i\ \p\<P* (1.4) 

where P T (E) denotes the probability of the event E with 
respect to the statistics SRB of the motions. 

Remarks: (1) In particular this holds for the single ob- 
servable a(x)/o~-\-, 0, EB|, and it is a theorem for hyper- 
bolic systems. This is the form in which the fluctuation 
relation was discovered, 0] , and the above is an extension 
of it under the same assumptions, (l3l. Il4j|. 

(2) The extension was found first in the special case 
/i = cr(x)/a + , / 2 = dE<y{x) where it has been shown 
to imply the Onsager reciprocity and Green-Kubo for- 
mulae, \u\ . 

(3) Another particular case is obtained by considering the 
probability that the averages f 1 , . . . , /„, over the time in- 
terval (— §j §) of the considered observables, have a given 
value ax, . . . , a n with p determined by a\, . . . , a n . Then 
for |p| < p* 

Pjfl ~ Ql> • ■ ■ ,fn ~ a n) _ e pa + r+0{l) g-j 
P(fl Cll, ---Jn a n) 

which is a very surprising relation because of the arbi- 
trariness of the observables fj which do not appear in 
the r.h.s. except through their function p. The above 
relation appeared recently in the context of Kraichnan's 
theory of passive scalars in a case in which p — X)i a «' 

M 

(4) A mahematically precise form of the theorem, 0, , 
is to say that for \p\ < p* the probability that p G A 
has the form exp(T max p6 A C(p) + 0(1)) and the func- 
tion C(p)i which for hyperbolic systems is known to be 
analytic and convex in a natural interval of definition 
(— p*,p*) (and —oo outside it) satisfies, for \p\ < p* the 
symmetry property 

a-?) = a?) - p<?+ (i-6) 
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(5) A further consequence is that the stationary average 
of exp J Q a{Stx)dt satisfies 

( e Io ^ StX '> dt ) ~i (1.7) 

where ~ 1 means that the quantity is bounded as r — ► oo 
(Bonetto's formula, [3, Eq.(9.10.4)]): see below for a 
hint to possible applications and for its similarity to 
Jarzynski's formula, [20l [2l|. 

(6) The chaotic hypothesis and the fluctuation relation 
have been tested quite extensively, starting with [f| , and 
it has almost become a test of the correctness of the 
computer programs simulations of chaotic systems rather 
than a formula to be tested. The situation is quite dif- 
ferent with experiments where a lot of difficulties arise, 
on a case by case basis, in setting up experiments and 
interpreting them. Nevertheless there have been sev- 
eral att emp ts and it can be hoped that more will come, 

(7) The formulae above hold for time evolutions described 
by maps iterations as well as for those described by dif- 
ferential equations. It is worth stressing that most works, 
in particular deal with discrete time evolutions. The 
case of continuous time is considered less frequently, and 
for the first time it has been formulated as a theorem 
in The case of maps is possibly closer to applica- 
tions as observations are usually done when some timing 
events occurr and evolution appears as a map between 
timing events. However a close examination of the rela- 
tions between the continuous and discrete cases reveals a 
number of delicate properties (particularly in the case in 
which singular forces may be acting, like Lennard-Jones 
type of repulsive cores) which if neglected may lead to 
errors, as exemplified in [28|. 

2. Entropy creation 

The result in the previous section hints at another ma- 
jor point of the research in the last 20 years of the '900's: 
in the '980's a concrete model of a thermostat became 
necessary to perform simulations of molecular dynamics 
in systems out of equilibrium. 

The "Nose-Hoover", the "isokinetic" or the "isoener- 
getic" thermostats are prominent examples that were put 
forward and employed to study a large number of prob- 
lems: see [29] for a thorough discussion of the related 
problems. One of the results was the discovery of the 
flutuation relation. 

Another important byproduct was the identification of 
the (not yet defined at the time) entropy creation rate 
with the phase space contraction. Although the original 
authors quite clearly attributed to the words they em- 
ployed a meaning close to the physical one suggested by 
the given names the general attitude was, it seems, to 
regard the thermostat models as unphysical and, conse- 
quently, to attribute little value to the concept of entropy 
creation as related in some way to the thermodynamic 
entropy. 



The above (extension) of the fluctuation theorem and 
relation, (1.4), suggests that one should give a, fundamen- 
tal physical sense to the phase space contraction, at least 
in finite, time reversible systems. Since as said above ulti- 
mately time reversal (or the equivalent, for our purposes, 
TCP symmetry) is a law of nature all models should ei- 
ther display the symmetry or be equivalent to symmetric 
models. 

In my view this identification between phase space con- 
traction and entropy creation rate is an important new 
development, [H, H3,[3l[, that is still not fully appreci- 
ated as it should. 

It has to be stressed that, although since more than 
a century we are familiar with the entropy of equilibria 
and its mechanical interpretation, no mechanical defini- 
tion of entropy creation rate in a process out of equi- 
librium has been proposed (or, better, accepted). The 
above extended fluctuation relation is clearly saying that 
the independence of the ratio in (1.5) means that the 
conditional probability that a pattern occurs in presence 
of an average (dimensionless) phase space contraction p 
is the same as that of the reverse pattern in presence of 
the opposite average phase space contraction (i.e. — p). 

In other words if the entropy creation rate is reversed 
during a time interval then evolution of the other observ- 
ables "proceeds backwards" with the same likelyhood it 
had to "proceed forward" when the average entropy pro- 
duction was p. All that has to be done to reverse the time 
arrow is the reverse the entropy creation rate. 

All this leads to say that the identification of entropy 
creation rate and phase space contraction has to be taken 
seriously. Its identity with what one would naturally call 
entropy creation in the many particular cases studied in 
the works summarized in [29J and continuing since should 
not be considered a curious coincidence but as a new 
insight into the foundations of nonequilibrium statistical 
mechanics. 

Note that we are saying that in nonequilibrium entropy 
creation is defined and identified with mechanical quan- 
tities: but entropy itself is not defined (yet?), not even in 
stationary states. It might even be not needed and not 
unambiguous, [HQ- 

To illustrate further the new identification of entropy 
creation rate with the mechanical quantity expressed by 
the divergence of the equations of motion, (1.3), and its 
physical interest I discuss a class of examples which, it 
seems to me, fully justify the mentioned identification, 

Consider a mechanical system Co in contact with the 
mechanical systems Cj, i — 1,2,.... Microscopically the 
state is is described by the positions Xq, Xi, . . . , X„ of 
the Ni, N±, . . . , N n constituent particles. The systems 
interact via short range pair forces with potential ener- 
gies J7 (Xj), C/ a (Xi,Xo). There is no direct interaction 
between the particles in C ai a > 0. 

The systems in C a , a > should be thought as "ther- 
mostats" acting on the system Cq across the separating 
walls via their mutual pair interactions. 
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^2 X , Xi, . . . , X, 







T X 0l = -a,(7o(Xo)-^] 9 j C/, s (X D ,X i )+F i 
X ai =-a i C/ Q (X a )-a i (7 a (Xo,X i )-aaX oi 

Fig.l: Schematic illustration of the geometry. The equations of 
motion are written here assuming unit mass for the particles. 

The thermostats temperature T a is defined to be pro- 
portional to the kinetic energy via the Boltzmann's con- 
stant k B . Setting K a = ±J2a X a *- f fa?T a N a , it is sup- 
posed constant and kept such by the action of suitable 
(phcnomenological) forces on the z-th particle in C a of 
the form — a a X a i. The a a can be taken 



ZNgk B Tg 



(2.1) 



where where the work performed by the system on the 
thermostats particles — J^X a • dx a U a (X.Q, X a ) can be 
called W a = Q a = heat given to the thermostat C a by 
the system in Cq. The external forces Fj are assumed to 
be purely positional. 

Given the above dynamical model (for heat trans- 
port) remark that it is reversible and time reversal is 
just the usual velocity reversal (because the thermostat 
forces are even under global velocities change). Fur- 
thermore the divergence of the total phase space vol- 
ume can be immediately computed and turns out to be 
cr(X) =cto(X) + /7(X) with 



to(X) 



E 

0=1 



Qa 3N a -l 



n 



Qa 



k B T a 3Ng 



^ k B T a 



(2.2) 



where U(X) = £"= 



3«„ 



and 0(N a 1 ) has been 



= 1 k B T a 3N a 

neglected in the last equality in (2.2). 

When computing time averages the "extra term" U 
will not contribute because being a time derivative its 
average will be ^([/(StX) — U(X.)) and therefore will 
give a vanishing contribution for large T and the average 

of J2a=i k^T = °o wm ^ e the average of a. If the inter- 
action U is bounded also the fluctuations of the averages 
of a and Co will coincide. In the cases in which the inter- 
actions are not bounded (e.g. Lennard- Jones repulsive 
cores) care has to be exercised in the fluctuations analy- 
sis: the picture does not change except in a rather well 
understood, trivial, way and this will not be discussed 
here, [Ull. 

Note that <tq is a "boundary term" , in the sense that 
it depends on the forces through the boundaries and the 
forces are supposed short range. Thus the question arises 
whether such kind of thermostats and short range inter- 
actions can lead to stationary states: this is not obvious 



but the "efficiency" of such thermostats has been investi- 
gated in molecular dynamics simulations, (35l.[36j, leading 
(not surpringly) to the result that the thermostat mecha- 
nism in Fig. 1 can lead to stationary states (even in pres- 
ence of additional positional forces stirring the particles 
in C ). 

Since the thermotats are regarded in equilibrium the 
above expression shows that <r(X) can be "legitimately" 
called the entropy increase of the reservoirs: so the me- 
chanical notion of phase space contraction acquires a 
clear physical meaning: and this is a no small achieve- 
ment of a long series of works based on simulations of 
molecular dynamics, (2^, [3l[ . 

3. Comments 

(1) Other studies of fluctuations have been proposed: 
they are rather different and apply to systems which are 
not stationary. The object of study are initial data sam- 
pled within an equilibrium distribution of a Hamiltonian 
system and subsequently evolved with the equations of 
motion of a dissipative time reversible system. 

Then the phase space contraction averaged over a time 

r, a d = i Jj" a(Stx)dt, will be such that the probability 
Po(a) with respect to the initial equilibrium distribution 
for a to have a given value is such that 



Po(a) 
Po(-a) 



(3.1) 



This is an exact identity, immediately following from the 
definitions. It involves no error terms, unlike the "simi- 
lar" (1.6) that can be written also as -^^y = e P a + T+ ° ( - 1 \ 
with P the probability with respect to the stationary dis- 
tribution, which is singular with respect to the equilib- 
rium distributions if <r+ > 0. 

It has been claimed that, being valid for all times, it im- 
plies the fluctuation relation, (1.6), for stationary states 
(at least when the stationary state exists). This would 
imply a simple, direct and assumptionless derivation of 
the fluctuation theorem in (1.7) and should hold in spite 
of the fact that in [4| an assumption about the chaotic 
nature of the motions is needed to derive it, together with 
a rather detailed understanding of the nature of chaotic 
systems. 

However a derivation of the fluctuation theorem (1.6) 
from (3.1) involves considering (3.1) after the limit r — > 
oo has been performed: a rather unclear procedure (note 
that the r.h.s. depends on r). Leaving aside the logical 
consitency problems it should be kept in mind that in the 
stationary state, at least in the interesting cases in which 
0+ > and there is dissipation, the statistics of motion 
will be controlled by a distribution that has nothing to 
do with the initial equilibrium distribution in which the 
averages in (3.1) are considered. 

Therefore the claim is incorrect and it is no surprise 
that some kind of chaos has to be present to obtain the 
fluctuation relation (1.6). In fact one can give examples 
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of simple systems in which (3.1) holds for all times, the 
system evolves towards a stationary state and neverthe- 
less the (1.6) does not hold, [37]]. The confusion has crept 
into the literature and even affected experiments: this 
can only be explained by a certain lack of attention to 
the literature due to the urge to find an easy way of test- 
ing the large fluctuations in real systems (fluids, granular 
materials, or even biological systems). 

Another aspect of (3.1) is that it involves a rather than 
p = a/a + . This is clearly a matter of convention: how- 
ever care has to be exercised because the fluctuation re- 
lation (1.6) is valid for \p\ < p* with p* > 1 being a 
physically nontrivial quantity, [19(. In terms of a this 
means that it is valid for \a\ < o~+. Overlooking this fact 
has led to think that it should hold for all a's and has 
led to errors in the literature. The errors are particularly 
noticeable in cases in which er + is close to zero, when 
the interpretation of simulations becomes quite difficult 
because long time scales become relevant. It has to be 
noted that cr^ 1 is a time scale diverging as cr + — * 0, see 
the discussion in (28j . 

(2) A different, interesting, fluctuation result is Jarzyn- 
ski 's formula which provides the means of computing the 
free energy difference between two equilibrium states at 
the same temperature. 

Imagine to extract samples x — (p, q) in phase space 
with a canonical probability distribution ^{dpdq) = 
ZQ 1 e~P H °( p > q ^dpdq, with Zq being the canonical partition 
function, and let So,t(p,q) be the solution of the Hamil- 
tonian time dependent equations p = —d q H(p,q,t),q = 

d p H(p^t)for < t < 1. Let H x (p, q) d = H{p, q, 1), 
then, [20, HH, 

Consider the "time 1 map" (p' , q') d = So,i(p, q) and call 

def 

W(p',q') = Hi (p , q ) — Hq (p, q) the corresponding vari- 
ation of the energy function. Then the distribution 
Z^ 1 e~^ Hlljp ' q ^dp'dq' is exactly equal to the distribution 
%i e -0W(p' V) fj, (dpdq). Hence 

, e -0W\ = ^1 = e -PAF W (3 2) 

where the average is with respect to the Gibbs distribution 
/io and AF is the free energy variation between the equi- 
librium states with Hamiltonians Hi and Hq respectively. 

Remark: (i) The reader will recognize in this exact iden- 
tity an instance of the Monte Carlo method. Its interest 
lies in the fact that it can be implemented without actu- 
ally knowing neither Hq nor Hi nor the protocol H(p, q, t). 
To evaluate the difference in free energy bewteen two 
equilibrium states at the same temperature of a system 
that one can construct in a laboratory, when the system 
changes its energy function from Hq to Hi (not necessar- 
ily explicitly known), then "all one has to do" is 
(a) To fix a protocol, i.e. a procedure, to transform the 
forces acting on the system along a well defined fixed once 



and for all path from the initial values to the final values 
in a fixed time interval (t — 1 in some units), and 

(b) Measure the energy variation W generated by the 
machines implementing the protocol. This is a really 
measurable quantity at least when W can be interpreted 
as the work done on the system, or related to it. 

(c) Then average of the exponential of —f3W with respect 
to a large number of repetition of the protocol abd apply 
(3.2). This can be useful even, and perhaps mainly, in 
biological experiments. 

(ii) Imagine a protocol consisting in lifting a container 
with a gas in equilibrium to height z: the Hamiltonian 
changes by Mgz, if M is the total mass and g gravity con- 
stant. Eq. (3.2) of course is correct, being an identity, 
and gives a free energy variation equal to f3Mgz while 
normally one would say that the free energy, and ev- 
ery other thermodynamic quantity, should have remained 
unchanged. Whether or not there has been a free en- 
ergy variation really depends on what one is interested 
in studying. Thus if the interest is in measuring free en- 
ergy variations in a biology experiment care has to be 
given (and is actually given) to make sure that the pro- 
tocol followed does not introduce spurious, quite hidden, 
forms of work. This makes, once more, clear that the 
application of a mathematical identity to real systems 
requires careful examination of the conclusions drawn. 

The two formulae (3.2) and (1.7) bear some similarities 
but are, however, quite different: 

(1) the f r a(S t x)dt in (1.7) is an entropy creation 
rather than the energy variation W. 

(2) the average in (1.7) is over the SRB distribution of 
a stationary state, in general out of equilibrium, rather 
than on a canonical equilibrium state. 

(3) the (1.7) says that (e ■>» B ( Sta dt ) SRB is bounded 
as t — ► oo rather than being 1 exactly unlike (3.2) which 
holds without corrections, [20l [2l|. 

The (3.2) has proved useful in various equilibrium 
problems (to evaluate the free energy variation when an 
equilibrium state with Hamiltonian Hq is compared to 
one with Hamiltonian Hi); hence it has some interest to 
investigate whether (1.7) can have some consequences. 

If a system is in a steady state and produces entropy 
at rate er+ (e.g. a living organism feeding on a back- 
ground) the fluctuation relation (1.6) and its consequence 
Bonetto's formula, (1-7), gives us informations on the 
fluctuations of entropy production, i.e. of heat produced, 
and (1.7) could be useful, for instance, to check that all 
relevant heat transfers have been properly taken into ac- 
count. This suggests that the fluctuation relation for 
stationary states could have some applications even in 
experiments in biology and be a valuable complement to 
(3.2). 

(3) Finally the identification of entropy creation and 
pase space contraction suggests a possible quantitative 
measure for "how irreversible" is a transformation be- 
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tween two different stationary states (equilibrium or not). 
Since physical processes are often accompanied by vol- 
ume changes with time V — > Vt it is natural to allow 
them and to change the definition of phase space con- 
traction (2.2) to, [331 ] . 

o x (K)=a (X) + U-N^ (3.3) 

Vt 

where N is the number of particles in the volume V. 

Then one can try to define the "irreversibility time 
scale" r(n) for a process II measuring the time scale 
over which the process manifests its irreversibility. Sup- 
pose that in the process II the parameters controlling the 
forces change with time from an initial value Fq to a final 
one Foe,; to be definite suppose 

F(t)=F + (l-e-T t )(F oo -F ) (3.4) 

The (3.4) allows us to consider the quasi instantaneous 
changes (7 — » 00) as well as the quasi static ones (7 — » 0). 

Starting the TV-particles system in the stationary state 
Ho with parameters Fo it evolves to fit and, eventually, 
to the stationary state (loo with parameters Fqo- 



Let Hsrb,t = be the SRB distribution with parameters 
F(t) "frozen" at value taken at time t. Then if af is 
the entropy creation rate in the "frozen" state fi S rb,t an 
irreversibility time scale for II could be defined as 

1 1"°° / \ 2 

r(Iirl = / (<°*V - <t**)sBB ,*) * (3-5) 

•J 

which can be checked to give the "expected results" in 
simple cases like the Joule expansion, see (33j . Quasi 
instantaneous processes II have a short t(T), mean- 
ing that irreversibility becomes noticeable immediately, 
while quasi static processes have a long r(n) indicating 
the opposite situation. 

(4) Although the above analysis is restricted to particle 
systems it can be extended to more general systems, in 
particular to fluids and turbulence. Not surprisingly as 
turbulence has been a source of inspiration for the devel- 
opment of the above ideas, fisl l33l|. 
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